Phase Transitions in Generalised Spin-Boson (Dicke) Models 
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We consider a class of generalised single mode Dicke Hamiltonians with arbitrary boson coupling 
in the pseudo-spin x-z plane. We find exact solutions in the thermodynamic, large-spin limit as a 
function of the coupling angle, which allows us to continuously move between the simple dephasing 
and the original Dicke Hamiltonians. Only in the latter case (orthogonal static and fluctuating 
couphngs), does the parity- symmetry induced quantum phase transition occur. 



I. INTRODUCTION 

Spin-Boson models appear in many areas of physics 
and are essential ingredients in theoretical quantum op- 
tics (light-matter interaction), nuclear physics 
quantum chaos 's^ , and quantum dissipation d] . The spin 
algebra can be used to describe single {j — N/2 —1/2) 
or many iV > 1 two-level systems where, in the sim- 
plest case, the interaction is with but a single bosonic 
mode (a, a^). Specific examples include cavity quan- 
tum electrodynamics and, more recently, 'phonon cav- 
ity quantum dynamics' of electrons interacting with sin- 
gle phonon (oscillation) modes in nano-electromcchanical 
systems 0, 0, S H such as freestanding quantum dots or 
'molecular transistors'. 

A common feature of spin-boson models is that in gen- 
eral they are non-integrable, with exact solutions avail- 
able only for very specific cases. Examples of the latter 
are simplified 'dephasing models', where the spin couples 
to both the boson and static field via only one of its com- 
ponents (usually chosen as Jz). Another example where 
exact solutions can be obtained is in the large spin limit 
00 where bosonic representations of spin Lie algebras 
[2| have been known for a long time; an early example 
being the Holstein-Primakoff transformation [3 • 

In this paper, we further explore the large-spin limit 
by starting from the most general, single- mode, spin bo- 
son Hamiltonian with linear coupling of all (a;, y, z) spin 
components to a static and a fluctuating (bosonic) term. 
For the specific case of the coupling of orthogonal [x and 
z) spin components to the static and the fluctuating term 
(Dicke model), we have previously found 'lOi IllL ll^ in- 
triguing connections between quantum chaos, entangle- 
ment, and the emergence of an instablity-induced quan- 
tum phase transition in the limit of large spins j ^ 00. 
Here, our main result will be that, surprisingly, this in- 
stability and the related parity-symmetry breaking of the 
ground-state wave functions only appears for 'orthogonal' 



coupling. The Dicke Hamiltonian [I^ (Rabi-Hamiltonian 
for spin 1/2) and its canonical equivalents therefore seem 
to be in a 'distinguished' class of Hamiltonians with very 
pronounced properties. It should be mentioned from the 
very beginning, however, that this distinction is most vis- 
ible in the strong coupling regime. 



II. THE MODEL AND ITS SOLUTION 



We start from a generic model Hamiltonian 



H 



(j7 + a^A + aAt) J, 



(1) 



describing the simplest coupling between Heisenberg- 
Weyl (l,a, flt) and the spin algebras Jx ~ + J-)^ 

Jy = i(J+ - J-), J2, with 



[J„J±]=±J±, [J+,J_]=2J, 



(2) 



In Eq. 1^ , O is a real and A a complex three-dimensional 
vector. Special cases of Eq.lpJ are the Rabi or the Dicke 
Hamiltonian [3 (rj = r^e^, A = At = Ae^), the 
simple dephasing Hamiltonian 0, 0, (fi = fie^, 
A = At = Aei with i = x,y, or z), the Jaynes-Cummings 
Hamiltonian 0(11 = flez, A = A[ex — isy]), and the 
one-mode version of the dissipativ e sp in-boson (tunnel- 
ing electron) Hamiltonian 0,11, EjEI = ^o^z+Te^, 
A = At = Ae^), where we denoted the unit vectors as 
Gi, i = x,y,z. The j = 1/2 variant of Eq. with 
Jl — wqGz + Tbx and A — ae^ + ibey appears in quasi 
one-dimensional quantum wires in the x-y plane in a con- 
stant magnetic field Bez for an electron gas with spin- 
orbit interactions (Rashba Hamiltonian) [20j. 

In the following, we restrict ourselves to A = At and 
therefore consider the Hamiltonian 



H = uja^a + ft ■ J + (a^ + a) A ■ 3 



(3) 



parameterised by two real three-dimensional vectors 
given by 
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(Aa;, Ay, Xz) 



(4) 



2 



where 1/^/27 is inserted to ensure correct scaling in the 
thermodynamic hmit, and the factor of 2 is for later con- 
venience. This Hamiltonian is invariant under a rotation 
about the z-axis, under which Jz and Jx — *■ —Jx, 

consequently, we shall only discuss the parameter range 
< < TT. Note that the more general case, Eq.(^, 
in general would leave three real, linearly independent 
three-dimensional vectors. The analysis would then be 
very similar to the following, though more cumbersome, 
which is why we restrict ourselves to the model Eq.Q. 

We begin by rotating our co-ordinate axes so that we 
work in the x-y plane, with the coupling-vector A aligned 
along the a:-axis. This gives us the form with which we 
shall work: 

2A 

H = uja)a + Vl{Jx cos + Jz sin 9) H — ■= (a^ + a) Jx (5) 
In deriving exact solutions for this model in the ther- 



modynamic limit, we shall follow the general procedure 
introduced for the Dicke model in Ref. [Tll |. 

First we employ the Holstein-Primakoff representa- 
tion of the angular momentum operators |9j, J-f = 
bW^j-b^b, J_ = v/2j - &t& b, Jz - (b^b-j). With 
Jx — ^ {J+ + J-), substitution gives us 

H = wflta + ^ cos6l (bW^j - 6+6 + ^/2j - Wb b^ 
+ 17 sin 6* (6+6 -j) (6) 

-I- — F= fa+ + a 

We next displace the oscillator modes a ^ a + ^Ja and 
b ^ b — where a and /? are assumed to be of the 
order of j. This leads to 



H = uj 

+ rising 



(a+a + V^{a^ + a) + a) + ifi cos oVk (fe+ 0? + 6 - 

(b^b - (6t + fe) + /3 - j) + A J A (at + a + 2^ (b^ ^ + Vv b ~ '^VPVv) , (7) 



r 



where 



6t6- V;3(6t + 6) 

fc = 2j-/3; r; = l ^} >- 

k 



(8) 



We now proceed to the thermodynamic limit, by taking 
j ^ oo and neglecting terms with powers of j in the 
denominator. This yields 



- ^ata+ Usin0 + 2AJ^ + ^^^y|j6t6+ Lv^_2AJ^ 



+ \4XJ^U- 13) ^nsme^ + ncose(^^)) {b^ + b) 



+n sin 9(p-j) + Lua-n cosOJkp ~ xj — (1 + ik) ~ - fl cos e\ -. 

y 2jk 4 V fc 



(9) 



r 



The two terms linear in bosonic operators can be elimi- and /3 is determined by 
nated by choosing the parameters a and /3 such that 



4A ak 



2A kP 
~V 2j' 



(10) 



^^ncoseVkfi-(pj^o. (11) 



3 



Substituting the value of a into this equation and sim- 
plifying, we obtain the following equation for ^/]3, 



4A2 J - (3 



J- [3 



0. (12) 



This equation is exactly soluble for y/j3, but the resulting 
form is extremely unwieldy. The solutions of this equa- 
tion for arbitrary parameters are cumbersome, and will 
not be reproduced here. In a few specific cases, to be elu- 
cidated later, compact expressions can be found. With 
the elimination of the linear terms, our Hamiltonian as- 
sumes the form 



H = uja'' a + Cjb'' b + s {b'' + b)'^ 

+ r (a^ + a) {b^ + b) + jEo + k' , 



(13) 



where the constants may be inferred by comparison with 
Eq. 0, with appropriate values of a and /3. Hamilto- 
nians of this form are analytically soluble via a unitary 
transformation, and since an example of this process was 
given in j^J, we shall not go into the details here. Suf- 
fice to say that after a Bogoliubov transformation of the 
bosonic operators, the Hamiltonian becomes diagonalised 



H 



+ jEG + k, (14) 

where we have introduced the excitation energies of the 
system, e±, and where Eq is the scaled ground-state en- 
ergy (scaled with j) and k is an unimportant constant of 
the order unity. In terms of the parameters introduced 
in Eq. H13() . the excitation energies are given by 

4 = (15) 
and, in terms of /3, the ground-state energy is given by 



JEg = n sine {(3 
2A_2 



+ 



3) 



(16) 



ncose^/f]{2j ~ (3). 



The general scheme in which we proceed from here is to 
solve Eq. 1)12(1 for (3, and then use this value to compute 
the excitation and ground-state energies. Before consid- 
ering the problem with arbitrary parameters however, we 
will focus on two special cases, which will explain many of 
the features of the general solution. It should be pointed 
out that not all solutions of Eq. H12|l are physically valid, 
and by considering the following cases we shall determine 
the criteria for selecting valid solutions. 



III. SPECIFIC LIMITS 

A. The Dicke model: 9 = tt/2 

In the case where the interaction and spin vectors are 
perpendicular we obtain the Dicke model: 

2A 

H^/2^^^a^a + ^Jz + ^^{a'' +a) J^. (17) 



In this limit there exists a conserved parity 11 such that 
[H, U] = 0, given by 



n = exp 



(18) 



where N is the "excitation number" and counts the total 
number of excitation quanta in the system. 11 possesses 
two eigenvalues, ±1, depending on whether the number 
of quanta is even or odd. 

For the Dicke Hamiltonian, the equation for determin- 
ing ^/P becomes 



(19) 



The simplest solution sets y/P — -Ja. = 0, which gives 
rise to the effective Hamiltonian 



H 



(1) 

72 



uiQb^b + uia^'a + A (a^ + a) (b^ + b) - jujo, (20) 



which has the excitation energies 



(cj2 - uj^y + WX^UJUJo } . (21) 



The excita- 



and the ground-state energy E'^' = ~joJa- 

tion energy e^-* remains real provided that A < Ac = 
and this demarcates the range of validity of 
this solution. The appearance of an imaginary part of 
an eigenenergy is one of our criteria for distinguishing 
between valid and invalid solutions of Eq. (|12|l . 

The remaining two solutions of Eq. H19|l are given by 
the displacements 



2A 



(l-M^), V^=±V7IW?, (22) 



where we have defined /x = = 4§- The Hamiltoni- 



K 

ans obtained with these solutions (one for each sign) are 
identical and have the same excitation energies 



n 2 



and ground-state energy. 



E. 



(2) 



2A2 



8A2 



+ 4Lu^Lun, (23) 



(24) 



and we thus see that these two solutions are completely 
degenerate. By considering the reality of e_ , we con- 
clude that these second two solutions are only valid pro- 
viding A > Ac. 

As described in |0 and to be discussed later, the exis- 
tence of these different solutions, one with zero displace- 
ment, and two with finite and opposite displacements, 
describes a quantum phase transition in the Dicke model, 
which occurs at the critical coupling Ac- The nature of 
this QPT is such that the parity symmetry becomes bro- 
ken above Ac, which explains the appearance of the two 
degenerate, broken symmetry solutions. 
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B. One-dimension: 6 — 

With interaction and spin aligned, the fuU Haniihonian 
of Eq. lO becomes 

2A 

+ 17 + —= ia^ +a)j^. (25) 

This Hamihonian is integrable for arbitrary j since its 
eigenstates are clearly also eigenstates of Jx, which al- 
lows us to replace the operator with its eigenvalue tUx — 
-j, -J + 1 • • • , j - 1, such that 



2A 

75! 



(a^ + a) ' 



(26) 



This leaves us with a single-mode bosonic Hamiltonian 
which may be diagonalised via a simple displacement a — )■ 
a — {2Xmx) I (V^^)- This results in the diagonal form 



Hq (j) — uja^a + Qnix 
which has the energy 



E„ 



flniT — 2- 



(27) 



(28) 



We proceed to the thermodynamic limit by writing rux = 
kx — j, and neglecting terms with j in the denominator. 
Whence, 



\ 2 \ / \ 2 - 

n + 4—]kx-jin + 2—\, (29) 

UJ J \ LO ' 



from which we immediately see that the excitation en- 
ergies are e_ = w and e+ = Vl + 4A^/w, and the scaled 
ground-state energy is Eq = — (p. -\- 2A^/w). 

We now seek to obtain this results using the general 
procedure outlined in section ^] The equation for the 
determination of /3 becomes 



(j - p) (p^uj + 4AV/3(2j-/3)) = 0. 



(30) 



Setting the second factor in this expression to zero leads 
to values of and \/a which give rise to complex exci- 
tation energies for all parameter values. These solutions 
are unphysical and we discard them as we did for the 
Dicke model. Considering the other solution, we have 
/3 = j, which gives y/P = iy^ and ^/a = ±^^/2j. With 
these choices, the Hamiltonian of Eq. (jl3|l becomes 



, ^\ At. 3 ^,2 1 



ijj 
2A^ 

UJ 



(31) 



Note that the two modes are now decoupled. The 6-mode 
may be diagonalised via the squeezing transformation, 



1 



(6t + ah) , 



1 



(5 + ah^) ,(32) 



with the squeezing parameter a = —1/3. In this way we 
arrive at the final form of the Hamiltonian 



The excitation energies of this Hamiltonian are clearly 
always real. However, only the Hamiltonian with the 
upper sign (corresponding to = has the same 

excitation and ground-state energies as our previous an- 
alytic calculation. The solution with ^/P = —^/j leads 
to a Hamiltonian with the incorrect energies, and is thus 
seen to be spurious. This solution is obviously unphysical 
for A < •\/wO/2, as here the coefficient of the second os- 
cillator becomes negative. The origin of this spurious so- 
lution can be easily understood by considering the 9 = tt 
limit of the Hamiltonian. In this case, the Hamiltonian 
is the same as that of Eq. H25I) , except that is replaced 
by — f2. Exchanging Jx for its eigenvalue as above and 
diagonalising the atomic mode, we obtain the energies 



E„ 



(34) 



The problem with this Hamiltonian arises when we take 
the thermodynamic limit under the assumption that 
nix — —j is the spin- quantum number of the ground 
state. This leads to the energy 

Kl^T =^n+ {-^ + 4AVt^) kx - j {-n + 2AV^) (35) 

which is the same as the spurious solutions obtained 
above. Clearly, the correct ground-state of the 9 — 
Hamiltonian actually has the quantum number nix — +J - 
So we see that the origin of this type of spurious solu- 
tion is due to the incorrect counting of the states labelled 
with nix as we go to the thermodynamic limit. The solu- 
tions with the incorrect sign always have a ground-state 
energy that is higher than the correct solution, and thus 
we are easily able to discard the solutions which arise 
from misidentifying the ground state. 



IV. RESULTS AND DISCUSSION 

To determine the behaviour of the system away from 
these two specific limits, we first solve for ^/a and \/]3. 
Figures ^ and |21 show the values of these two displace- 
ment parameters as functions of both A and 9. Our first 
observation is that ioi all 9 ^ ■k/2, there is only one solu- 
tion for a given A. Furthermore, the sign of y/a and 
is given by that of cos 9. The divide between the regions 
of positive and negative displacements is spanned by the 
special case of = 7r/2, which is the previously discussed 
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FIG. 1: The two displacement parameters ^Ja and \/^ as a 
function of the couphng A for various different angles 6. The 
Hamiltonian is on scaled resonance, = f2 = 1, Ac = 0.5 



-- X=0.5 
i=0.2 




FIG. 2: The two displacement parameters ^/a and \f]5 a func- 
tion of the Q for representative couplings. The Hamiltonian is 
on scaled resonance, uj = Q, = \, \c = 0.5 

Dicke model. In this case we have ^Ja — — below 
Ac, and two solutions of opposite sign above Ac. The dis- 
placement parameters a and /3 determine the centre(s) of 
the collective ground-state wave function of the coupled 
systems in a position-momentum representation of the 
two bosonic modes a and h [lol |. The appearance of two 
solutions for = 7r/2 then corresponds to a breaking up 
of the wave function into two macroscopically separated 
parts for j — > cxd. This parity breaking phase transition 
therefore occurs only at = tt/2 which demonstrates 
that the Dicke model with its 'orthogonal' coupling 
is unique within the whole class of Hamiltonians Hg. It 
is only in this special case that the super-radiant phase 
will exhibit macroscopically coherent (Schrodinger's cat) 
behaviour when j remains finite. 

This conclusion is corroborated by considering the ex- 
cited states of our models. The nature of the system is 
characterised by the behaviour of its two excitation en- 
ergies, which are plotted in Figs. |31and21 In Fig. |2|the 
limiting cases of = and 6* = 7r/2 are clearly identifi- 
able, and serve to provide bounds for the other solutions 




FIG. 3: The excitation energies of the system as a function of 
the coupling A for various different angles 9. The Hamiltonian 
is on scaled resonance, lj = Q — 1, \c = 0.5 
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FIG. 4: The excitation energies of the system as a function of 
angle 6 for representative values of coupling A. The Hamilto- 
nian is on scaled resonance, u> = = 1, Xc — 0.5 

away from these values. The most crucial consequence of 
this is that again, only for 9 = 7r/2 and A = Ac does e_ 
identically vanish, and so it is only for these parameter 
values that a quantum phase transition occurs. 

A further check is made in Fig. [Sj where we plot the 
values of important observables of the system. The ex- 
pression for the ground-state energy has been given in 
Eq. (|24|) . The atomic inversion and mean field occupa- 
tion are given by 

{Jz)/j^P/j-i; {a^a)=a/j. (36) 

Again, singular behaviour in the form of non-analyticities 
of the curves at A = Ac = 1/2 is observed only at 6* = tt/2 
in agreement with the above result. 

To summarise, the existence of the quantum phase 
transition for spin-boson models Hg is dependent on the 
two vectors A and being exactly perpendicular, which 
one might not have expected at the outset. In conclusion, 
we briefly discuss the implications these findings have for 
spin-boson systems. One obvious consequence is that 



6 




'non-orthogonal' coupling terms always would smear out 
phase transitions or their precursors when tuning from 
a weak to a strong coupling regime in, e.g., photon or 
phonon cavities. At first sight, this looks like bad news 
for the possible realisation of critical behavior in real- 
istic systems where one would always expect perturba- 
tive terms leading to a general, not necessarily orthogo- 
nal coupling, unless some symmetry prevents this from 
occuring. On the other hand, it would be desirable to 
explore tunable systems where one can vary the parame- 
ter 9 (for example by using external electric or magnetic 
fields), in order to test some of our predictions. 



FIG. 5: The ground-state energy, atomic inversion and mean 
photon number of the ground state as a function of the cou- 
pling A for various different angles 6. The Hamiltonian is on 
scaled resonance, tJ = f2 = l, Ac = 0.5 
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